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EMPIRICAL MODELING OF THE THERMOSPHERE: AN OVERVIEW 
A. Hedin, NASA/Goddard Space Flight Center 


Hedin gave a summary of thermospheric density modeling 
history and standard atmospheres. In particular, he compared 
and contrasted the approaches of the Jacchia and MSIS models. 

His conclusions were that the Jacchia models are best if drag is 
the primary concern. MSIS is superior for variations in 
composition and temperature variations and comparison with 
theoretical models is facilitated by the use of spherical 
harmonics, which also provide a simple and consistent way of 
obtaining simplifications. 

ADVANTAGES / DISADVANTAGES 

1. Jacchia 

a. Theoretically best if drag is primary quantity 
desired without high resolution and for satellite 
geometries and orbits similar to those used in 
generating the model. However, drag coefficients 
used in density derivation need to be more 
carefully specified if original drag is to be 
reproduced. Inaccurate specification of compo- 
sition (e.g. He bulge) may result in inaccurate 
drag. 

b. Absolute total density dependent on the drag 
coefficient rather than the instrument 
calibration. However, dependence of drag 
coefficients on composition and extreme 
geometries may be a problem. Model predictions of 
composition and temperature are derived from 
auxiliary data or assumptions and may not be 
realistic. 

c. Formulation has particular difficulty coping with 
minor constituent variations found by mass 
spectrometers and cumbersome pseudotemperatures 
of J77 help only a little. 

2 . MSIS 

a. Best for composition/temperature variations, but 
agrees with drag models in overall averages. 

b. Provides better resolution of variations 
(including total density) in local time, etc. 
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Absolute densities dependent on individual 
calibration constants for contributing instruments 
but model accuracy should be better than that of 
an individual instrument. 

Spherical harmonics facilitate systematic 
increase in model resolution and comparison with 
theoretical models. Similarly, complexity can be 
reduced if desired by dropping higher harmonics 
or unneeded effects. 
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MODEL DEVELOPMENT 


A. Jacchia 

1. J65 

a. Earliest comprehensive model based on drag. Lower boundary at 120 
loo. 

b. Static height profiles a3 function of exospheric temperature 
assuming hvrfro<»*-.»M ^/diffusive pnn 1 1 < brliim 

c. First to include four principal effects (diurnal/seasonal, 
semiannual, solar activity, magnetic activity) using ad hoc 
formulas for exospheric temperature to fill gaps. 

d. Introduced Bates type temperature profile (which can be integrated 
explicitly) . 

2. J70 i J71 

a. Lower boundary at 90 km and more complicated temperature profile 
requiring numerical integration. 

b. Refinements and expansions of ad hoc formulas. 

c. Included factor of three winter helium bulge. 

d. J71 raised atomic oxygen at 150 km over J70. 

3. J77 

a. Inclusion of some results from mass spectrometers. 

b. Magnetic coordinates for magnetic activity effects. 

c. Composition phase through pseudo- temperatures. 

B. OGO-6/MSIS 

1. OGO-6 (1974) 

a. Earliest comprehensive model based on mass spectrometer data. 

b. Bates temperature profile above 120 km. 

c. Spherical harmonics for geographical/local time coordinates. 

d. Variable boundary at 120 km for He and 0 to represent phase 
differences between constituents. Height profiles assuming 
hydrostatic/diffusive equilibrium. 

e. Temperature inferred from N2 agreed well with incoherent scatter. 

2. MSIS 77 

a. Same format as 0G0-6. 

b. Used mass spectrometer density data from five satellites and 
temperatures from four incoherent scatter stations. 

c. Variable boundary also for N2 so temperature depends on incoherent 
scatter and N2 scale heights. 

3. MSIS 79 

a. Introduced UT/Longitude variations for quiet and magnetic active 
times (alternative to magnetic coordinates). Temperature maximum 
and He minimum near magnetic pole. 

4. MSIS 83 

a. Density and temperature data from mass spectrometers on seven 
satellites, from five IS stations, and from rockets. 

b. Extended profiles below 120 km to 85 km using analytically 
integrable temperature profiles. 

c. Includes major variations in temperature and density below 120 km. 

d. Improved resolution in prediction of magnetic activity variations 
using time history of 3hr indiciea. 
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CONCENTRATION (cm 



LATITUDE (°) 

Fig. 13 

Latitudinal variation of n {He) at 1 000 km altitude. The left 
part corresponds to F = F — 150 x 10 -22 W m“ 2 Hz -1 
and Kp = 2. The right part corresponds to F = F = 92 x 

10 -22 W m - ^ Hz” 1 and K = 2. Comparison with 
Jacchia 1971 and MSIS models. Barlier (1979) 
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EXOSPHERIC TEMPERATURE DISTRIBUTION AT NORTHERN SUMMER SOLSTICE 

Jacchia J71 Modal 
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Figure 8. Exospheric isotherms (°K) above the globe, computed from equa- 
tions (15) and (16), for the case when the minimum temperature 
is 1000°K. Aitoff's equal-area projection; meridians of local solar 
time and parallels of latitude are drawn 30® apart. Top, equinoxe 
bottom, northern summer solstice. 
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Table 4a. Density Ratio to ISIS-83 for N^, 0, and He 




h 



0 



He 


Data Set 

altitude 

avg sd 

pts 

avg 

sd 

pts 

avg 

sd 

pts 

OGO-6 (HS) 

400-700 

1.08 .27 

659 

1.15 

.16 

1276 

1.18 

.19 

902 

San Marco-3 (IS) 

190-250 

1.10 .20 

77 

.86 

.15 

24 

1.09 

.17 

41 

Aaroa-A HATE (MS) 

200-500 

1.13 .47 

321 

1.14 

.33 

478 

1.18 

.42 

466 

AE-C NATE (MS) 

190-400 

1.13 .33 

640 

.91 

.18 

866 

.68 

.18 

855 

AE-C OSS (MS) 

135-1604- 

.97 .15 

440 







AE-C OSS (MS) 

190-400 

1.02 .26 

319 

1.08 

.18 

387 

1.03 

.23 

371 

AE-0 OSS (MS) 

140-1604- 

.99 .16 

184 







AE-D OSS (MS) 

190-400 

.87 .33 

99 

1.01 

.18 

107 

.78 

.22 

107 

AE-E NACE (MS) 

140-1604- 

1.01 .13 

815 







AE-E NACE (IS) 

190-450 

1.00 .22 

701 

.87 

.18 

1019 

.93 

.17 

1002 

ESRO-4 (MS) 

200-350 

.88 .33 

427 

.83 

.24 

587 

.84 

.30 

518 

Rockets (MS) 

100-120 

.83 .36 

35 







Rockets (MS) 

110-160 

.92 .30 

28 







Rockets (MS) 

190-300 

.90 .32 

39 







Areclbo (IS) 

100-120 

.92 .32 

228 







Arecibo (IS) 

110-135 

1.14 .51 

109 
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